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We determine accurate values of ordering temperatures and critical exponents for Ising Spin Glass 
transitions in dimension 4, using a combination of finite size scaling and non-equilibrium scaling 
techniques. We find that the exponents r] and z vary with the form of the interaction distribution, 
indicating non-universality at Ising spin glass transitions. These results confirm conclusions drawn 
from numerical data for dimension 3. 



I. INTRODUCTION 

Precise values of critical exponents at Ising Spin Glass 
(ISG) transitions have been hard to obtain from numer- 
ical simulations principally because of the agonizingly 
slow relaxation near the transition. Finite size scaling |^] 
can help bypass this problem but can introduce other 
difficulties, associated with the possible need to allow 
for corrections to the finite size scaling rules. We have 
shown that in favourable cases a combination of finite 
size scaling and non-equilibrium scaling techniques can 
allow one to estimate to high accuracy both the ordering 
temperature and the exponents 77 and z. We exploited 
this method for ISGs in dimension 3 where we found that 
the exponents r/ and z changed considerably as a func- 
tion of the distribution of near neighbour interactions, in 
complete contradiction to what would be expected from 
standard second order universality rules. There is no 
formal theorem which states that ISG transitions should 
obey universality rules, but up to now general sentiment 
tended to be in favour of universality. The empirical ob- 
servation of non-universality Q implies that the standard 
application of the renormalization group method to ISGs 
must be reconsidered, and shows that there are impor- 
tant qualitative differences between ISG transitions and 
conventional second order transitions. 

We have now applied the same approach to ISGs in di- 
mension 4 and for one case in dimension 5. We find that 
in dimension 4 as in dimension 3 the exponents rj and z 
vary systematically from system to system. The numer- 
ical results for the exponents as functions of dimension 
can be compared with the predictions of the e expan- 
sion; agreement is very poor if the expansion is extended 
beyond the leading term. 

II. METHOD 

Technically the problem of measuring exponents is less 
arduous in higher dimensions than in dimension 3. Be- 
cause dimension 4 is well above the lower critical dimen- 
sion the existence of a finite Tg in dimension 4 has never 



been in doubt, in contrast to the situation in dimension 
3. The Binder cumulant curves for each system show 
clear intersections and there already exist quite precise 
estimates of the ordering temperatures Tg, at least for 
the ±J and Gaussian interactions 

We will first outline the simulation techniques that we 
have used. First we measure, for a number of test tem- 
peratures T near the probable value of Tg, the time de- 
pendent non-equilibrium spin glass susceptibility X*(i). 
Two completely random (infinite temperature) replicas 
A and B of the same system are quenched to T with 
independent updating. The spin glass susceptibility 

X*it) = [< SHt)Sf{t) >2] (1) 

is recorded as a function of time t after the start of the 
quench. Precisely at Tg, x* increases with t as t^, with 
/i = (2 — 1^. If we do not know Tg a priori, we obtain 
an apparent h(T) at each test temperature T . This non- 
equilibrium parameter presents the considerable advan- 
tage of requiring no preliminary anneal, and represents 
the growth of correlations as the internal temperature 
drops in the ISG samples. Because we concentrate on 
short times after having chosen an appropriately sized 
sample the correlation lengths are always much smaller 
than the sample size so there are no finite size correc- 
tions to a very good approximation. This type of non- 
equilibrium scaling has been tested very carefully on reg- 
ular systems where the ordering temperature and the ex- 
ponents were already known |^ , and the scaling has been 
rigorously verified. 

Secondly, at the same test temperatures T we anneal 
samples for a waiting time tw before measuring the initial 
decay of the autocorrelation function over a further time t 

q{t)^<S,{t^)S,{t + t^)> (2) 

If t « and T = Tg, q{t) = with x ^ {d - 
2 + r])/2z Under these conditions the decay of q{t) 
is identical to the equilibrium form of the decay (which 
would be measured after an infinitely long anneal time) , 
and again if the sample size is chosen appropriately there 
are no finite size corrections lOl. It turns out that that 
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at temperatures T < Tg, q{t) follows an algebraic decay 
law with effective values of x which are smaller than the 
value corresponding to T = Tg. For higher temperatures 
there is a further multiplicative factor /{t/r). 

It can be seen that if reasonable precautions are taken 
the effective exponent combinations h{T) and x{T) can 
be measured with negligable systematic error. One can 
thus obtain a first set of effective exponents 771 (T) and 
z{T) as functions of the test temperature T from the 
equations 
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g the normalized spin glass susceptibility xsg/L 
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[0 giving an independent measure 0.10 



Finally we use finite size scaling of the equilibrium spin 
glass susceptibility as a function of sample size L. Again 
at T, 

proportional to L^^ 
of 77. Below Tg the power law form continues to hold, 
providing a second measure of 77, 772 (T). Now if we plot 
?7i (T) and 772 (T) against T, consistency dictates that the 
true Tg and 77 must correspond to the intersection point 
of 7,1 (T) and 772 (T). 



made on samples of size L = 20 with anneals of 10^ MCS 
and runs of Iff^ MCS. Runs were taken with one sample. 
Data points for h{T) and x{T) are shown in figure ||. 
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FIG. 1. The effective h{T) and x{T) (see text) as functions 
of temperature T for the ±J interaction (o), the Gaussian 
interaction (□) and the decreasing exponential interactions 
{<>). Full symbols: x. Empty symbols: h. 



III. RESULTS IN DIMENSION 4 

We have made measurements for ISGs on the 4 dimen- 
sion simple cubic lattice, with near neighbour interac- 
tions of different types. In practice the finite size scaling 
is the most demanding part of the simulation as far as 
CPU hours are concerned, as it is essential to make strin- 
gent checks that equilibrium has been achieved, and av- 
erages must be made over a large number of samples. We 
have been fortunate in being able to use equilibrium XSG 
data from A. P. Young |3] for the 4d ±J interaction case 
and from Parisi et al |5| for the 4d Gaussian interaction 
case. For the other two cases we carried out simulations 
in the usual manner for L up to 8 taking standard pre- 
cautions to be certain of complete thermal equilibrium 
at each size. 

We thus dispose of data for the 4 dimension ±J (J), 
Uniform (U), Gaussian (G), and decreasing exponential 
(Ed) interactions. Earlier work had shown that Tg 
is close to 2.0 for the ±J case and close to 1.75 for the 
Gaussian case. To get a trial estimates of Tg for the Ed 
and U cases we have relied on a Migdal-Kadanoff ap- 
proach where it is assumed that the normalization 
scaling parameter b is the same for all members of a fam- 
ily of ISGs in a given dimension. This gives quite reliable 
values of Tg once the value for one member of the family 
is known. Simulations for h{T) were made using up to 
2000 samples of size L = 10. x{T) measurements were 



It can be seen that for each system the effective value 
x{T) increases with T while the effective value h{T) de- 
creases with T. In dimension 3, h{T) increases with T 
we have no explanation for this difference. 
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FIG. 2. The effective dynamic exponent z obtained using 
equation ^ for ±J interaction, gaussian interaction and de- 
creasing exponential interactions, (same symbols as figure |l|). 
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We first extract z{T) from the h{T) and x{T) data 
using equation |3[ figure ||. Clearly z is not universal; 
z tends to increase steadily from one type of system to 
the next as the kurtosis of the interaction distribution 
increases and Tg drops. 

In figure ^ we have plotted 771 (T) and r]2 (T) for the ± J 
case where we have the most complete r]2 data. 
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FIG. 3. Effective exponents ??i(T) (o) and 772 (T) (•) (see 
text) as function of temperature for the ± J interaction case. 
The intersection gives Tg and the exponent rj. 

There is a clear intersection at the point Tg — 2.00 ± 
0.01 and r; = -0.30 ± 0.02. High quahty Binder cu- 
mulant data for L = 4 to 12 due to Young |^ show 
a clearly defined crossing of all the gL{T) curves at 
T = 2. 00 ±0.01 |l4|. Two series expansion calculations 
gave Tg estimates of 2.02 ± 0.06 § and 2.04 ± 0.05 §. 
The agreement between different determinations of Tg in 
this case is thus very satisfactory, giving further confi- 
dence in the new method outlined above. We can note 
that for estimating 77, this method is fairly insensitive to 
the exact value of Tg because the line for rji (T) is much 
less steep in the region of Tg than is the line for 772 (T) . 
Also the present method estimates 77 at the ordering tem- 
perature only, and so does not imply any assumption of a 
scaling relation over a range of temperatures around Tg. 

If we now turn to the other sets of interaction distri- 
butions. Figure ^, we find similar intersections leading to 
the Tg and rj values given in Table |. 

The Tg value in the Gaussian interaction case lies be- 
tween the two published estimates obtained from Binder 
cumulant work, and is close to estimates we can obtain 
starting from the ± J value for Tg and using the Migdal- 
Kadanoff technique [T^ or by using the series expansion 
formula due to Singh and Fisher |15| (both methods give 
essentially Tg 



to be 2.00 for the ±J case) 



1.2 



1.0 



0.8 



0.6 



0.4 



0.2 





1.40 



1.50 



1.60 



1.70 



1.80 



1.90 



T 



FIG. 4. As in figure |^, for the Gaussian interaction case (□) 
and the decreasing exponential case (O). Empty symbols: 771. 
Full symbol: 772 

For the U and Ed cases there are no Binder cumu- 
lant data to compare with but the Migdal-Kadanoff and 
Singh-Fisher techniques lead to Tg values of 1.90 and 1.53 
respectively for the two interaction distributions. Over- 
all agreement between the simulation estimates and these 
values can again be seen to be excellent. The 77 values. 
Table |, like the z values, are strikingly non-universal, 
with 77 becoming more negative as the kurtosis of the 
distribution increases and the Tg drops. 



IV. DIMENSION 5 

For the ±J interaction ISG in dimension 5, series ex- 
pansion gives a precise value for the ordering tem- 
perature, Tg = 2.57 ± 0.01. We have seen above that 
in dimension 4 the simulations and the series work led 
to very similar Tg estimates, so the dimension 5 series 
estimate should be very reliable. We have asumed this 
Tg value is correct and have measured h and x at this 
Tg. The values. Table " 



lead to z = 4.50 ±0.1 and 
?7 = —0.39 ± 0.02, in excellent agreement with the series 
estimate 77 = -0.38 ± 0.07. 

We have also measured h for the Gaussian ISG in di- 
mension 5 {Tg = 2.31), and for the ±J ISG in dimension 6 
at Tg which is equal to 3.03 |@,|6l. 



1.75 for the Gaussian case if Tg is taken 
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V. CRITICAL EXPONENTS - DIMENSIONAL 
DEPENDENCE. 

We can first concentrate on the behaviour of the crit- 
ical exponents and their combinations as a function of 
dimension for the series of ISGs with ± J interactions, in 
order to compare with e expansion expressions Jl7| , p^ . 
In the standard renormalisation group treatment of the 
ISGs, the upper critical dimension is 6 and calculations 
have been made to third order in e {e ~ 6 — d). 

1/ 77(d). The values of 77 from Q and the present work 
are shown in Figure ^ with an extrapolation to the upper 
critical dimension value 77 = for dimension 6. 
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FIG. 5. ?7 as a function of dimension for the ± J interaction 
(o), the Gaussian interaction (□) and the decreasing exponen- 
tial interactions (O). The full lines are to guide the eye. The 
straight dashed line indicates the e expansion to first order; 
the dot dashed curve is the e expansion to order 3. 

The renormalisation group e expansion estimate to 
leading order in e and to order 3 is also indicated. It can 
be seen that the simulation values, which from the dis- 
cussion given above we consider very reliable, behave in 
a regular fashion. The initial trend of 77(0?) extrapolated 
towards d = 6 is broadly consistent with the e expansion 
to lowest order. However the e expansion curve to order 
3 lies a long way from the numerical points. This is in 
striking contrast to the standard second order transition 
case where the e expansion to the same order gives ex- 
cellent agreement with numerical or analytic values right 
down to e = 2. 

2/ h{d). Using van Hove arguments, Zippelius p^ ] 
found that that below the upper critical dimension z = 
2(2 — 7;), with no correction to leading order in e. From 
the definition of h, this relation can be simply rewritten 
h = 0.5, so we would expect h to be close to this value 



at and near to d = 6. Indeed we find that h is equal to 
0.5 at d = 6, and that h(d) does stay close to 0.5 down 
to dimension 4, Figure ^. 
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FIG. 6. The exponent ft as a function of dimension for the 
±J interaction (o), the Gaussian interaction (□) and the de- 
creasing exponential interactions (O). The point at dimension 
6 is measured. 

This shows that the Zippelius relation between z and 
77 is fairly accurate. 

3/ z(d). At dimension 6 we expect 2 = 4; the z data as 
a function of dimension are approaching 4 as dimension 
is increased towards 6, Figure ^. 
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FIG. 7. The dynamic exponent z as a function of dimension 
for the ±J interaction (o), the Gaussian interaction (□) and 
the decreasing exponential interactions {O). The point at 
dimension 6 is the theoretical upper critical dimension value. 
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4/ h'{d). We have no new information on v from our 
simulations, but for completeness we present results com- 
piled from various sources PJ^JTo|l in Figure ^ together 
with the e expansion curve. 



are 2 — d and infinity respectively. If we plot, Figure ||, 
ry as a function of Tg for dimensions 3 an 4, we find that 
the trend of the numerical values is consistent with a 
tendency for 77 towards 2 — d as Tg tends towards zero. 
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FIG. 8. A collection of critical exponents for the ±J inter- 
action case as a function of dimension. The v values are taken 
from references p|j7|jic|] and the 7 and fJ values are calculated 
from 1] and u using the standard scaling relations. The full 
points at d = 6 correspond to the theoretical upper critical di- 
mension values and the dashed lines represent the theoretical 
e expansions to first order. 

Once 1] and v are known, the other static exponents 
can all be deduced using the scaling rules. Figure ||. 
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FIG. 9. The trend of ?7 as a function of Tg for dimension 3 
(circles) and dimension 4 (squares). 

Plotting, Figure |l^, 1/z against Tg we again find a 
trend corresponding to a divergence of z as Tg tends to 
zero. 



0.25 



VI. CRITICAL EXPONENTS - INTERACTION 
DEPENDENCE. 

It can be seen in Table | that (as in dimension 3), the 
exponents in dimension 4 vary in a systematic manner, 
with Tj becoming more negative and z tending to a higher 
value when the kurtosis of the interaction distribution in- 
creases and the Tg drops. (There are many other possible 
parameters which could be used to modify the interaction 
distribution - for instance an alternative way to increase 
the kurtosis would be to dilute the interactions. Wc do 
not know if there is a one to one relationship between the 
exponents and the kurtosis in all cases). 

We can attempt to understand from a nai've stand- 
point what the systematic behaviour is indicating [ p^ . 
If we make the heuristic assumption that the kurtosis 
is a pertinent parameter as far as the critical exponents 
are concerned, then as the kurtosis is increased and the 
Tg is driven towards zero, we would expect 77 and z to 
tend concomitantly towards their zero Tg values, which 
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FIG. 10. i plotted against Tg for dimension 3 (□) and 
dimension 4 (o). i extrapolates toward zero as Tg decreases, 
indicating that z diverges when Tg goes to zero. 

We have so far studied the effect of just one parameter, 
the form of the interaction distribution, on the critical 
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exponents in finite dimension ISGs. Other parameters 
readily suggest themselves; it may well be that there is 
a continuous variation of the values of the critical expo- 
nents whenever any parameter is changed, so that the 
universality class concept simply does not apply to spin 
glasses. It is certainly relevant that the critical exponent 
z has been shown analytically to change continuously as 
a function of applied field along the AT line for the mean 
field spin glass . It appears that the physics of tran- 
sitions in complex systems such as spin glasses is funda- 
mentally different from that of second order transitions 
in regular systems. 

VII. CONCLUSION 

We have found accurate values for the the ordering 
temperatures Tg and for the critical exponents r\ and z 
in Ising Spin Glasses with different sets of near neighbour 
interactions on the simple cubic lattice in dimension 4. 
The results demonstrate that, as in dimension 3 ||^, the 
exponents are not universal, confirming that the univer- 
sality class concept is not relevant to spin glass transi- 
tions. This observation shows that the spin glass transi- 
tion cannot be treated as a second order transition in the 
conventional sense, and that the renormalization group 
theory as applied to spin glass transitions (and by ex- 
tension to transitions in all complex systems) should be 
seriously reconsidered. 
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d System Tg h{Tg) -q z__ 

4 ±J 1.99±.01 0.51±.02 -0.31±.01 4.45±.l 
U 1.91±.01 0.51±.02 -0.37±.02 4.65±.l 
G 1.77±.01 0.48±.02 -0.47±.02 5.10±.l 
Ed 1.52±.01 0.43±.02 -0.60±.03 6.05±.l 

5 ±J 2^57 0.530±.02 -0.39±.02 4.50±.l 
G 2.31 0.537±.02 — — 



TABLE I. Values of the temperature of transition and crit- 
ical exponents for ISG for various distributions of interactions 
in dimension 4 and 5. The temperature of transition in di- 
mension 5 where taken from series expansion 
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